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ABSTRACT 

We generalize the spherical collapse model for the formation of bound objects to apply 
in a Universe with arbitrary positive cosmological constant. We calculate the critical 
condition for collapse of an overdense region and give exact values of the characteristic 
densities and redshifts of its evolution. 

Key words: methods: analytical - cosmology: theory - galaxies: clusters general - 
galaxies: formation - large-scale structure of Universe 



1 INTRODUCTION 

The spherical collapse model was first developed by Gunn 
& Gott (1972) for a flat Universe with no cosmological con- 
stant. It assumes that the process of formation of bound 
objects in the Universe can be at first approximation de- 
scribed by evolution of an uniformly overdense spherical re- 
gion in otherwise smooth background (and it is therefore 
called the top hat model). Despite its simplicity, the model 
is still widely used to explain properties of a single bound ob- 
ject via extensions such as the spherical infall model (Gunn 
1977; Hoffman & Shaham 1985; Lokas 2000) as well as sta- 
tistical properties of different classes of objects via Press- 
Schechter-like formalisms (Press & Schechter 1974, hereafter 
PS; Lacey & Cole 1993, 1994; Sheth & Tormen 1999). 

Recently, our knowledge on background cosmology has 
improved dramatically mainly due to new supernovae and 
cosmic microwave background data. Current observations 
favor an almost fiat Universe with Qq = 0.3 (see e.g. Harun- 
or-Rashid & Roos 2001 and references therein) and the re- 
maining contribution in the form of cosmological constant 
or some other form of dark energy. A considerable effort has 
gone into attempts to put constraints on models with dark 
energy (Wang et al. 2000; Huterer & Turner 2001). Another 
direction of investigations is into the physical basis for the 
existence of such component with the oldest attempts going 
back to Ratra & Peebles (1988). One of the promising gen- 
eralizations of the cosmological constant is the quintessence 
model (Caldwell, Dave, & Steinhardt 1998) based on so- 
called "tracker fields" that display an attractor-like behavior 
causing the dark energy density to follow the radiation den- 
sity in the radiation dominated era but dominate over mat- 
ter density after matter-radiation equality (Zlatev, Wang, & 
Steinhardt 1999; Steinhardt, Wang, & Zlatev 1999). 

Given the growing popularity of models with cosmo- 



logical constant we generalize the description of the spher- 
ical collapse to include its effect both in flat and non-flat 
cosmologies. The top hat model serves as a basic tool in 
performing analytic calculations of structure formation via 
gravitational instability in an expanding Universe, most no- 
tably in the framework of the PS formalism. Our aim here 
is to extend the arsenal of analytical, or quasi-analytical, 
formulae describing the redshifts and (over) densities char- 
acterizing the coUapse processes to the case of a Universe 
dominated by a cosmological constant. We derive some sim- 
ple analytical formulae and flts that will serve as useful tools 
in constructing models of structure and galaxy formation. 

The paper is organized as follows. In Section 2 we briefly 
summarize the properties of the cosmological model with 
cosmological constant including the linear growth factor of 
density fluctuations. Section 3 is devoted to the evolution of 
the overdense region and gives the critical threshold for col- 
lapse. Sections 4 and 5 discuss the characteristic densities of 
the forming object and redshifts of evolution. The summary 
follows in Section 6. 



2 THE COSMOLOGICAL MODEL 

The evolution of the scale factor a = R/Ro = 1/(1-1-2:) (nor- 
malized to unity at present) in a Universe with cosmological 
constant is governed by the Friedmann equation 

da _ Ho 
dt " 7W 
where 

-1/2 



(1) 



/(«) = 



l-F!^o(i-l) +\o{a-l) 



(2) 



and Ho is the present value of the Hubble parameter. 
The quantities with subscript here and below denote the 
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present values. The parameter is the standard measure 
of the amount of matter in units of critical density and A 
measures the density of cosmological constant in the same 
units 

^= — =3^ (3) 

where A = const is the cosmological constant in standard 
notation. 

The evolution of Q and A with redshift z is given by 

2 

(4) 



n{z) = no(i + z) 

and 

X{z) = Ao 
where 



Ho 
H{z) 



Ho 
H(z) 



(5) 



H{z) 



Ho 



= no{l + zf - (Qo + Ao - 1)(1 + zf + Ao. (6) 



The linear evolution of the matter density contrast S 
Sg/g is governed by equation 



S + 2-5-4TTGgS = 



(7) 



where dots represent derivatives with respect to time. In the 
presence of cosmological constant the growing mode can be 
constructed in a simple form (Heath 1977; Carroll, Press, & 
Turner 1992) 



D{a) = 



5^0 



f{a)da 



(8) 



'2a f (a, 

where /(a) was defined in equation (^). The expression in 
(^) is normalized so that for Q = 1 and A = we have 
D(a) = a. 

For some special cases one can obtain analytical expres- 
sions for D{a). In the well-studied case of Ao = and Qo < 1 
we have 



D{a) = 



a 



3/2 



2F1 



- - 1 ) S^o + 1 

a 



1/2 



(9) 



3 5 7 
2' 2' 2' 



which is equivalent to the better known expressions given by 
e.g. Peebles (1980) or Padmanabhan (1993). For Qo + Xo = 1 
we get 



a 



3/2 



2F1 



An ^ ■ 
I a I i2o -I- a 

5 3 11 3/, ^\ 
6' 2' 6 V no) 



(10) 



For arbitrary (rio,Ao) pairs D{a) is easily obtained by nu- 
merical integration in equation (H) (see also Hamilton 2001). 



3 EVOLUTION OF THE OVERDENSE 
REGION 

We assume that at some time ti corresponding to redshift 
Zi the region of proper radius ri is overdense by the aver- 
age Ai = const with respect to the background, that is it 
encloses a mass 



M(n) = yPb,irf(l + AO (11) 

where ph.i is the background density of matter at ti. 

Evolution of this region is governed by the energy equa- 
tion 

'dr\'^ GM H^Xr^ _ 



1 

2 \dt 



GM 
r 



Expressing e as a combination of the kinetic and potential 
energy per unit mass at ti, 

e{tO = ^[l-m + ^0-Xi], (13) 

using equation ( ]ll| ) and introducing a new variable s = r/ri, 
equation (|l^) can be rewritten in the form 

ds _ Hi 
dt g{a,s) 

where 



(14) 



1+^1(1 + A; 



+ Ai(s'- 1) 



-1/2 



(15) 
A(2;i) are 



The parameters Hi = H{z\), fii — Q.(zi), Ai 
given by equations (^)-(^). 

Since the energy in the overdense region is conserved 
we can use equations (p^-(p^) to determine the maximum 
expansion (or turn-around) radius rta (or equivalently, Sta = 
rta/ri) of the overdense region. It must obey the following 
condition 



felSta + ^'2Sta + h-j, = Q 

where 

fei = Ai 

62 = l-f^i(l + Ai)-Ai 

63 = a(l-hAi). 

There is one interesting solution to equation ( p^ 

2 /-h2\^^'^ /(I>-2ti- 



Sta — 



^/3 V 61 / V 3 y 



where 
S — arccos 



(a;2 2/2)1/2 
with 

X = -961/^63 



2ni1/2 



y = [3(-4fe^ - 276i6^)] 

For Ao = we simply get sta = — 63/62. 

The condition for the solution ( |l8| ) to exist is 

Ai > Ai,c. = j^^(Ai) - 1 
where 

s , 5Ai 3Ai(8 + Ai) 3v{\i) 
It Ai = 1 + -r + ' ^ + — 
4 4t;(Ai) 4 

and 

V = v{Xi) = {Ai[8 -\- + 20Ai + 8(1 - Ai)^''^]}'/^ 



(16) 
(17) 

(18) 

(19) 

(20) 
(21) 

(22) 
(23) 
(24) 



In the limit of Ao ^ we have u{\i) 1. In this limit we 
reproduce the well known condition for the overdense region 
to turn around (see e.g. Padmanabhan 1993) 
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Figure 1. The ratio Aj cr/^i cr(-^0 = 0) as a function of Qq 
for Zi = 100. Solid curves from top to bottom show results for 
Ao = 0.2,0.4,0.6,0.8 and 1. The thin horizontal line marks the 
limiting case of Aq = 0, where the two overdensities are equal. 



T 1 1 1 1 1 1 1 r 




Figure 2. Evolution of radius of the presently collapsing over- 
dense region in models with and without cosmological constant. 



Ai > Ai,„(Ao = 0) = — - 1. 



(25) 



The condition equivalent to ( P4 ) has been discussed 
previously by Weinberg (1987), Martel (1991) and Lokas 
& Hoffman (2001). For all positive Ao, we have u(Xi) > 1 
so the condition (E2) seems more stringent than (Eq), i.e. in 
the Universe with cosmological constant overdensities have 
to be larger in order to collapse than in the Universe with the 
same density parameter Q,i but no cosmological constant. It 
is more useful, however, to compare the conditions for cos- 
mological models with the same value of the present density 
parameter, f^o- Since the evolution of O depends on A (see 
equation and Q), the relation between conditions ( |22[ ) 
and ( psj ) for a given ilo is not obvious. It turns out that Ai,cr 
can in fact be higher as well as lower than Ai,cr(Ao = 0) de- 
pending on the choice of Qo and Aq. Figure |i| shows the ratio 
of the two critical overdensities calculated for 21 = 100 as a 
function of Qo for different values of Ao — const. 

It is worth noting that in the case of Ao = condition 
( psi ) is equivalent to the requirement of the energy of the 
overdense region to be negative. In the case of a Universe 
with positive Ao, according to equation (^^ the condition 
_E < translates into Ai > ^^"^(l — Ai) — 1 which allows 
densities much lower than (123) since (1 — Ai) < 1. 

Integrating equation ([13) numerically for given fio and 
Ao we get the trajectory r{t). Figure]^ compares examples 
of r{t) obtained for f^o = 0.3 and two values of Ao. For the 
Qo = 0.3, Ao = case we have the well known analytical 
solutions r/rta = (1 — cos6')/2 and i/icoii = {0 — sin 60/(271) 
with Q < 9 < 2-K. The plots of r{t) shown in Figure H were 
obtained with the assumption of tcoii equal to the present 
age of Universe. In the Ao = case r{t) is independent of 
the collapse time, while calculations for Ao 7^ show that 
r{t) approaches the Ao = solution for ^coii 00. 



4 THE CHARACTERISTIC DENSITIES 

4.1 The linear density contrast at collapse 

Integrating equations (^) and ( |l4[ ) we get 



f{a)da = Hot 



g{s)ds = Hit. 



(26) 
(27) 



Eliminating t we obtain equations which can be used to 
calculate the scale factor at turn-around (ata) and collapse 
(fflcoii) of a region with particular Ai at a given zi 



f(a)da = ^ / 9{s)ds 
Jo 

"coll rr /"sta 

f{a)da = 2-^ / gis)ds 







Hi 



(28) 
(29) 







where Sta is given by equation ( |l^ ) and we defined the col- 
lapse time tcou to be twice the turn-around time 



ita = y f{a)da. 



(30) 



Assuming that the mass inside the overdense region 
does not change, the overdensity inside the sphere of size 
r with respect to the background density at any time is 



Pb s-^ 



where we used equation (pT|). 



(31) 
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At early times, t 0, we can expand the expressions 
on the left-hand sides of equations (|2^)-(^^ around a = 
and s = respectively. Integrating term by term we obtain 



Hot = 
Hit = 
+ 



3/2 , no + Ao-1 5/2 , ^,7/2 



3/2 



3[ni(l + Ai)]l/2' 



50, 

.^3/2 



a"'-' +0{a''') (32) 
(33) 



ni(l + Ai)+Ai-l 5/2 ^.7/2, 

5[f7i(l + Ai)p/2 ^^^^^ 



Inverting both series we express a and s as power series of 

^2/3 



,2/3 I ,4/3 I ^/,8/3\ 
a = at ' + C2t ' + U(t ' ) 

where 

/■?\2/3 

c, = (^) (i^o^no)^/^ 

3 /3\i/' ^^<y'(l-f^o-Ao) 



(34) 



C2 



1/3 



and 

S = dlt"'" +d2t'''-' +0{t 

where 

,3s 2/3 



(35) 



(I) 



3_ 1^3^/3 Hf/=^[l-»i(l + Ai)-Ai] 
10 



[a(l + Ai)]i/3 

Inserting expressions ( ^^ and (^Hj) into equation (|3l| 
and keeping only the lowest order term we find 



3 /3\2/3 



H^^''{l^no-Xo) 



2/3 



(36) 



//f''^[fii(l + Ai) + Ai-l] 



+ 



[f7i(l + Ai)]2/3 



4/3n 



From expansion ( p^ ) we have to the lowest order t^'^^ 
so the dependence of 5 on a in the limit of a is 

S = h{Qo, Ao, Ai, 2i)a + O(a^), 

where 



/i(r2o, Ao, Ai, Zi 
_ 3 

^ 5 



(37) 
(38) 



1 - fio - Ao _^ [Q.i{l + Ai) + Ai - 1](1 + zi) 



no 



ni(l + Ai)2/3 



Using equations (y)-(|6|) we can rewrite h in terms of the 
present cosmological parameters 



/i(f2o, Ao, Ai,Zi 

_ 3 

^ 5 



(39) 



I -no- Xo _^ no[l + Ai(l + Zi)] + Ao - 1 



^0 ' no(l + Ai)2/3 

For arbitrary (Slo, Ao), one can show that 
D{a) = a + 0{a^) 



(40) 



is the properly normalized solution of equation (|7|) around 
a = 0. Given this behavior of the linear growth factor D(a), 
we finally obtain the density contrast as predicted by linear 
theory 
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Figure 3. Upper panel: parameter <5c as a function of f2o 
with solid lines corresponding from bottom to top to Aq = 
0, 0.2, 0.4, 0.6, 0.8, 1 and the dashed line showing results for the 
flat case + Aq = 1. Lower panel: <5c as a function of Zj-oii 
for four models with Ho = 0.3 and Ao = 0,0.4,0.7,1. The thin 
straight line marks the fiducial value 5c = 1.68647 for f2o = 1, 
Ao = 0. 



5l = /i(!:!o, Ao, Ai,2:i)_D(a). 



(41) 



Note that D{a) is the general solution to equation ([^. The 
first order expansion of D{a) has been used only to derive 
/i(no,Ao, Ai,Zi). 

A particularly useful quantity is the linear density con- 
trast at the moment of collapse i.e. when s reaches zero 



5c = /i[f2o, Ao, Ai(acoii), zi]i?(acoii). 



(42) 



Ai(acoii) in the above equation means that Ai corresponding 
to acoii has to be determined for a given Zi from equation 
(p^. The problem can be reduced to solving equation ( |29| ) 
with Sta given by equation (^sj). 

Figure 1^ shows 5c for different values of Aq. In the upper 
panel the assumption is that the collapse is taking place 
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now, Zcoii ~ 0. Solid lines display the dependence of 5c on 
Qq with Ao = const while the dashed line has Ao chosen so 
that r2o + Ao = 1. In this last special case we reproduce the 
results of Eke, Cole & Frenk (1996), while for open models 
with no cosmological constant our results match those of 
Lacey & Cole (1993). 

The lower panel of Figure ^ shows how 5c changes with 
the redshift of collapse, Zcoii, for a few models with Qo = 0.3 
and different values of Aq. We see that the dependence on 
ZcoU is rather weak and at large Zcoii the values converge to 
the well known fiducial value of 5c = 3(1271)^/^20 ^ 1.68647 
valid in the Universe with fio = 1 and Ao = 0. This value is 
particularly quickly reached with increasing ZcoU for the flat 
case, flo = 0.3, Aq = 0.7. 

4.2 The density of virialized halo 

Another useful quantity is the ratio of the density in the 
object to the critical density at virialization 

f2(acoii) 



pvir 



(flc 



(^) [l + Ai(aeoii)] (43) 



Pcrit 

where Scoii = '"coii/^'i and rcoii is the effective final radius of 
the collapsed object. We assume that the object virializes 
at tcoii, the time corresponding to s ^ 0. Application of the 
virial theorem in the presence of cosmological constant leads 
to the following equation for the ratio of the final radius of 
the object to its turn-around radius F = rcou/rta. (Lahav et 
al. 1991) 



where 



(2 + 77)F+l = 



2AisL 



(44) 



(45) 



a(i + Ai) 

In the calculations of Ac we use the exact solution to equa- 
tion (^) which in the case of A > can be written down 
using expression ( |l^ ) with F instead of Sta and bi = 27], 
b2 = —(2 + 77) and — 1. However, a good approximation 
is provided bv F « (1 - r?/2)/(2 - ri/2) (Lahav et al. 1991). 

Figure H shows Ac for different Aq. The upper panel 
gives its values as a function of Qo for models with different 
values of Ao with the assumption that the collapse occurs 
at ZcoU = 0. The solid lines correspond to constant values of 
Ac while the dashed line has Aq chosen so that f^o + Aq = 1. 
Again, we agree with the results for the special cases of 
fio < 1, Ao = and fio -I- Ao = 1 derived previously by 
Lacey & Cole (1993) and Eke et al. (1996) respectively. 

We see that for a given Qo lower Ao makes virialized ob- 
jects denser with respect to critical density. It is interesting 
to note, however, that this relation is inverted for higher ZcoU 
as proved by the lower panel of Figure ^ where we display 
the dependence of Ac on the redshift of collapse for a few 
models. Again, as in the case of 5c, we observe that at high 
^coii values of Ac approach the well known fiducial value of 
Ac = IStt^ « 177.653 valid for fio = 1, Ao = and the 
convergence is fastest for the fiat case. 



5 CHARACTERISTIC REDSHIFTS 

It is sometimes useful to be able to estimate the redshift 
of a particular stage of evolution of the perturbation given 
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Figure 4. Upper panel: parameter Ac as a function of f2o 
with solid lines corresponding from top to bottom to Ao = 
0,0.2,0.4,0.6,0.8,1 and the dashed line showing results for the 
flat case Ho -|- Ao = 1. Lower panel: Ac as a function of Zcoii 
for four models with Ho = 0.3 and Ao = 0,0.4,0.7,1. The thin 
straight line marks the fiducial value Ac = 177.653 for Ho = 1, 
Ao = 0. 



its present overdensity as predicted by linear theory, 5l (a = 
1) = (5o. For the redshift of collapse combining equations 
(|l|) and (H) we obtain 



5q = 5c(aca\\) 



D{a ^ 1) 
-D(acoii) 



(46) 



Using the previously obtained results for 5c and appropriate 
formulae for the linear growth of fiuctuations (^), (|l^) or 
(^), we can calculate the present linear density contrast of 
fiuctuation that collapsed at Zcoii. This relation can only be 
inverted analytically in the case of fio = 1, Aq = when we 
get Zcoii = 5o/5c — 1. For other cases the calculations have 
to be done numerically. 

Using equation analogous to ( ^^ we can also calcu- 
late the redshift of turn-around, zta.. Parameter 5c(acoii) has 
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Table 1 . Values of the best-fitting parameters a and /3 of equation 
@ for Ho = 0.3, Ao = 0.7. 



2 


a 


/9 


accuracy 


Zcoll 


0.774 


1.15 


3 % 


2ta 


1.25 


1.30 


5 % 


•2nl 


2.32 


1.15 


3 % 



then to be replaced by the corresponding turn-around value 
5ta(ata) which we do not give here, but which is calculated 
numerically from equation (|^. 5ta is close to unity for all 
models considered here. 

Another interesting epoch in the evolution of an over- 
dense region is the onset of nonlinearity, which we character- 
ize here by redshift Zni- This is the time when the nonlinear 
density contrast given by equation ( |3l| ) reaches unity. Again 
equation equivalent to maj can be used with Sc replaced by 
5ni(ani), which turns out to be of the order of 0.5. We obtain 
Sni from equation ( p^ replacing ata and Sta by a^i and s^i 
where Oni and Sni obey equation (|l|) with S — 1. 

Figure ^ presents the characteristic redshifts Zcoii, 
and Zni as functions of So for different models with Qo = 0.3 
and cosmological constant Ao — 0,0.4,0.7 and 1. The thin 
solid line in each panel gives for reference the exactly linear 
relation in the fiducial case of f2o = 1, Ao = 0. 

Results shown in Figure ^ suggest that the dependence 
of characteristic redshifts on So can be well fitted with simple 
linear formulae 



z — a So ^ (3 



(47) 



where different constants a and /3 correspond to each of the 
three characteristic redshifts. In the f2o = 1, Ao = model 
this relation is exact, we have /3 = 1 and Qcoii ~ l/5c, Qta = 
l/(5ta, Qni = l/^ni- For the most popular flat Universe with 
Slo = 0.3 and Ao = 0.7 we find the best fitting parameters a 
and (3 shown in Table ^. The fits are intended to be useful at 
high redshifts where the relation between z and So is almost 
exactly linear. The accuracy of the fits in terms of z obtained 
for a given So for z > 1 is shown in the last column of the 
Table. It should be emphasized that for high redshifts the 
fitted values of a and (3 work much better than Sc, (5ta, Sni 
and 13 = 1. 



6 SUMMARY 

The top hat model has been extended here to the cases of 
arbitrary positive cosmological constant including non-fiat 
cosmologies. In particular, we have calculated the critical 
(over)density for collapse, the virial density and the charac- 
teristic redshifts of the collapse process. These include the 
redshifts of the transition to nonlinearity, turn-around and 
the collapse epoch. The characteristic redshifts cannot be 
represented by closed analytical expressions and therefore 
simple fitting formulae have been provided. 

The top hat model constitutes the basic tool used in an- 
alytical and semi-analytical models of large scale structure 
and galaxy formation. The prime example here is the calcu- 
lation of the mass function of collapsed objects by the PS 
formalism and its recent extensions (Lacey & Cole 1993; 
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Figure 5. Redshift of collapse (upper panel), turn-around (mid- 
dle panel) and transition to nonlinearity (lower panel) of density 
fluctuation with present linear density contrast 5o for four models 
with Co = 0.3 and Aq = 0,0.4,0.7 and 1. The thin straight line 
in each panel marks the fiducial case of Qo = 1, Ao = 0. 
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Somerville & Kolatt 1999). The PS mass function has a 
Gaussian term of the critical (over)density for collapse and 
therefore an exact evaluation of this density should be used. 
Another important application of the top hat model is the 
semi-analytical modelling of galaxy formation (Kauffmann 
et al. 1999; Somerville and Primack 1999). The virial pa- 
rameters of collapsed objects and the properties of baryons 
within such objects are often calculated in the framework of 
the top hat model, and this depends crucially on the back- 
ground cosmology. 



Evolution of spherical perturbation 7 



ACKNOWLEDGEMENTS 

This work was supported by the Polish KBN grant 
2P03D02319 and by the Israel Science Foundation grant 
103/98. 



REFERENCES 

Caldwell R. R., Dave R., Steinhardt P. J., 1998, Phys. Rev. Lett., 
80, 1582 

Carroll S. M., Press W. H., Turner E. L., 1992, ARAA, 30, 499 

Eke V. R., Cole S., Frenk C. S., 1996, MNRAS, 282, 263 

Gunn J. E., 1977, ApJ, 218, 592 

Gunn J. E., Gott J. R., 1972, ApJ, 176, 1 

Hamilton A. J. S., 2001, MNRAS, 322, 419 

Harun-or-Rashid S. M., Roos M., 2001, A&A, 373, 369 

Heath D. J., 1977, MNRAS, 179, 351 

Hoffman Y., Shaham J., 1985, ApJ, 297, 16 

Huterer D., Turner M. S., 2001, Phys. Rev. D, 64, 123527 

Kauffmann G., Colberg J. M., Diafcrio A., White S. D. M., 1999, 

MNRAS, 307, 529 
Lacey C, Cole S., 1993, MNRAS, 262, 627 
Lacey C, Cole S., 1994, MNRAS, 271, 676 

Lahav O., Lilje P. B., Primack J. R., Rees M. J., 1991, MNRAS, 
251, 128 

Lokas E. L., 2000, MNRAS, 311, 423 

Lokas E. L., Hoffman Y., 2001, in Spooner N. J. C, Kudryavtsev 
v., eds, Proc. 3rd International Workshop, The Identifi cation 



of Dark Ma tter. World Scientific, Singapore, p. 121, astro- 
ph/( |011295| 



Martel H., 1991, ApJ, 377, 7 
Padmanabhan T., 1993, Structure formation in the Universe. 

Cambridge Univ. Press, Cambridge 
Peebles P. J. E., 1980, The Large-Scale Structure of the Universe. 

Princeton Univ. Press, Princeton, NJ 
Press W. H., Schechter P., 1974, ApJ, 187, 425 (PS) 
Ratra B., Peebles P. J. E., 1988, Phys. Rev. D, 37, 3406 
Sheth R., Tormen G., 1999, MNRAS, 308, 119 
Somerville R. S., Kolatt T. S., 1999, MNRAS, 305, 1 
Somerville R. S., Primack J. R., 1999, MNRAS, 310, 1087 
Steinhardt P. J., Wang L., Zlatev I., 1999, Phys. Rev. D, 591, 270 
Wang L., Caldwell R. R., Ostriker J. P., Steinhardt P. J., 2000, 

ApJ, 530, 17 
Weinberg S., 1987, Phys., Rev. Lett., 59, 2607 
Zlatev I., Wang L., Steinhardt P. J., 1999, Phys. Rev. Lett., 82, 

896 



© 0000 RAS, MNRAS 000, 000-000 



